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Abstract. A generalized theory of multi-field galileons has been recently put forward. This 
model stems from the ongoing effort to embed generic galileon theories within brane con- 
structions. Such an approach has proved very useful in connecting interesting and essential 
features of these theories with geometric properties of the branes embedding. We investigate 
the cosmological implications of a very restrictive multi-field galileon theory whose leading 
interaction is solely quartic in the scalar field ir and lends itself nicely to an interesting cos- 
mology. The trispectrum of curvature fluctuations has features which are quite distinctive 
with respect to their P(X, (p) counterpart. 

We also show that, despite an absent cubic Lagrangian in the full theory, non-Gaussianities 
in this model cannot produce the combination a of small bispectrum alongside with a large 
trispectrum. We further expand on this point to draw a lesson on what having a symmetry 
in the full background independent theory entails at the level of fluctuations and vice- versa. 
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1 Introduction 

Quantum field theories endowed with galileon symmetry are ubiquitous in recent literature. 
They were first seen at work in the context of brane constructions, specifically as a limit 
of the DGP [1] model and have been investigated ever since. One can study their higher 
dimensional origin or be content with an effective 4D theory, obtained integrating out the 
bulk, which inherits most interesting properties. It is this last route that we will follow here 
but we refer the interested reader to the literature for the former perspective: [1-3] . 

Among the most compelling properties which originate from the galilean symmetry 
characterizing these theories are second order equations of motion (from a theory with higher 
order derivatives) and non-renormalization theorems that protect the coefficients of galileon 
terms from large renormalization. The first property translates into a theory with a well 
defined Cauchy problem which is possibly unitary in the quantum regime and the second 
into a model which is predictive as, under suitable conditions, the number of terms which 
describe the theory stays finite and their coefficients are only subject to small corrections. 

It would then be intriguing to see if any such theory can accommodate for an inflationary 
mechanism. Recently, an interesting concrete realization of such a model, galileon inflation, 
has been proposed in [6]: it represents one of the very few radiatively stable and properly 
predictive models of inflation to date. This model shares with DBI [4] inflation essential 
properties such as being unitary and having symmetries at one's disposal which protect 
leading operators from large quantum corrections. 

DBI inflation provides us with a very compelling inflationary mechanism embedded in 
a UV finite theory such as string theory. However, at least for its single- field-disguise, DBI 
predictions are quite close to be ruled out by data; the Planck mission is likely to narrow 
the bounds on cosmological observables which already look unfavorable for this model (not 
so for multi-field models, see e.g. [5]). 

The results of [6] are compatible with current observational constraints: the bispectrum 
of curvature perturbation has been studied and predictions for the observable Jnl have been 
extracted. The calculation of the corresponding trispectrum has been performed in [16] (see 
[17] for a different approach). As per usual, given a model for primordial perturbations, one 
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would like to obtain the corresponding predictions on all available observables, chief among 
which are non-Gaussinities (NG) [28]. Ideally, NG of a given model are within experimental 
values and present distinctive features in the form of the so-called higher order correlators 
amplitudes and shape- functions. 

Considering the overwhelming amount of inflationary models which are in agreement 
with observations and predict more or less the same scenario for CMB observables, one, while 
waiting for further data, might want to look for some guiding principles in the search for the 
most compelling inflationary theories. 

Below we elect the stability of the theory as paramount and introduce a stable inflation- 
ary model which has quite distinctive signatures at the level of the trispectrum of curvature 
fluctuations. Our model is based on the co-dimension 2 (quite in general, it would suffice 
to have co-dimension 2n) galileon theories introduced in [3] and the way we implement an 
inflationary mechanism is in close correspondence to what is done in [6]. 

The authors of [3] show how an even co-dimension forces an additional symmetry 
(SO(N)) on the galileon theory which propagates down to the 4D effective action and results 
in a Lagrangian with a quadratic and quartic piece only. As it turns out, when perturbing 
around a de Sitter background, the third order action for perturbations is not parametrically 
suppressed, this despite a null £3 in the full theory. A small-bispectrum — large-trispectrum 
scenario cannot therefore be implemented. This prompts some comments on how to imple- 
ment such a scenario in the full theory. 

The paper is organized as follows: in Section 2 we review the model first introduced in 
[3]. In Section 3 we show how such a model can implement an inflationary phase. Section 4 is 
devoted to the calculation of non-Gaussianities. We offer some comments in the Discussion 
section, then summarize our findings and possible future work in the Conclusions. 

2 Review of the Model 

Here we describe in some detail the model first discussed in [3] . The key feature that identifies 
this model is the additional (with respect to, say, [2]) SO(N) symmetry it endows to the 
galileon fields. This symmetry emerges from the requirement that the theory behaves as its 
co-dimension 1 counterpart and it is therefore both connected with the essential interesting 
properties of galileon theories and extremely natural. 

Following [3], we will take below a relatively quick route to the galilean action of in- 
terest; for a more systematic higher dimensional derivation of the same theory, we refer the 
reader to [3]. 

In flat space, upon requiring that a theory is 

• second order in the equation of motion and 

• endowed with galilean symmetry ir — > it + c + b^x^, 

one arrives [29] to the following general formula: 

C n+ i = T^yy ^(-1)V 1 ^V 2P( ^ ) -..V 1 " pK) (d fll ird Ul vd tl2 d v ^...d» n d Un K) ; (2.1) 
V p 

there are d non trivial terms (plus the tadpole) in d space-time dimensions, five in 4D. Sys- 
tematics ways of deriving the corresponding equation of motion for tt and currents associated 
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with the shift symmetry have been already explored in the literature, and we do not report 
the results here. What concerns us is the generalization of Eq. (2.1) to the multifield case: 



Cn+1 = (n 2 -l)\ E(- 1 )V 1P( "V 2PM ---^ PK) {^ n+1 d^d Vl ^d^d V2 K I \..d tln d v y™) , 
v / ■ p 

(2.2) 

where 5 , / 1 / 2 .. is a tensor symmetric in the field index /. It is easily verified that this multifield 
generalization preserves the galilean symmetry on each field tt 1 — > tt 1 + c + b^x^ and retains 
equations of motion free of the Ostrogradski [9] instability. As we will see, both in the single 
and multifield case the galileon symmetry can be though of as deriving from symmetries 
which characterize a higher dimensional brane construction [2]. 

Consider indeed the case of a 3-brane in a higher dimensional bulk (5D for now). Assume 
flat bulk and use X A (x) to describe the embedding of the brane, where A signals the bulk 
dimensionality and x describes the brane coordinates. We want any action to be invariant 
under Poincare transformations of the bulk and gauge invariant under reparametrization of 
the brane: 

S P X A = co^X B + e A ; S g X A = ^(x)d^X A (2.3) 

where e a describes translations in the bulk, ui Lorentz transformations and £ is the usual 
gauge parameter. One can fix a gauge, e.g. unitary gauge, at this stage: 

X»{x) = x», X b {x) = tt(x) => 5 P X» = u»x v + u^vr + e", (2.4) 

and so it becomes clear that an additional input is needed to have a gauge fixed action 
invariant under the Poincare transformation. Upon choosing: 

= -uj£x u - w£tt - e", (2.5) 

the combined action 5 pi = 5p + 5 g transformation is a symmetry the gauge fixed action. The 
crucial bit now is to see the effect of this transformation on ir which, the notation gives it 
away, is going to be our galileon. Indeed the total transformation on tt reads [6]: 



5 p/7 r = -u^x^tt - e^vr + - u^d^ + e 5 (2.6) 

The first two terms are the unbroken 4D Poincare transformations, the second two terms are 
the broken boosts and the last one corresponds to broken translations in the fifth direction; 
in other words what happens is I SO (1,4) —> I SO (1,3) . To get a glimpse of the galilean 
symmetry consider the internal relativistic invariance under which tt: 

5 pi7 r = - w^vt^tt + e, (2.7) 

where the index 5 has been omitted. In the non relativistic limit (small 7r) this is the usual 
tt — > ir + c + b^x^, but we now see how this symmetry is coming from the brane motion in 
the (flat) bulk. This results lends itself to a straightforward generalization. Suppose that 
the co-dimension is not 1 anymore, but generic. The only relevant difference for us will be 
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on the transformation for tt, now tt in the gauge fixed action under the combined effect of 
5 pm 1 = dp-K 1 + SgTr 1 gives: 

5pm 1 = -^ u x u dy - e^d^ir 1 + wjc" - w> J d M 7r 7 + e 1 + JjTt 3 . (2.8) 

Let us focus on the last term, it is the only one which is not a straightforward generaliza- 
tion of Eq. (2.6). Indeed uj^ is antisymmetric and so its o;| entry vanishes in the single field 
case. This last term is signalling us an SO{N) symmetry in the dimensions traverse to the 
brane, therefore on top of the usual galilean invariance on each field, in the non relativistic 
limit we now count an additional SO{N). 

This symmetry has striking consequences on the galilean Lagrangian of the type in 
Eq. (2.2). Indeed, with Si x j a ,„ a symmetric tensor, as a matter of simple indices contractions, 
one realizes the only allowed actions contain an even number of fields, in 4D we have just 
two different contributions! These are the usual kinetic term and a quartic term in the field 
tt: 

C 2 = d^d^TT 1 , (2.9) 
U = d^dy (d^d p 7T J d u d p 7T J - d^d u 7T J nir J ) + ^d^d^ir 1 (□vr J n^ J - d v d p % J d v d p TT J ) , 

where the indices /, J are contracted with a Kronecker delta. Having reached this point, 
it is useful to remark that both the renormalization and (second order) equation of motion 
properties of this multifield theory (the latter can be seen basically by inspection) are left 
untouched. It is the galilean symmetry that protects these theories from large renormalization 
corrections and such a symmetry is still intact. 

The additional SO(N) symmetry is the news and, as we will see in next section, it can 
be put to good use in cosmological scenarios. We note here that perhaps the most elegant 
route to the Lagrangian above is the one that arrives at it by writing down the most generic 
action with galilean symmetry and internal relativistic invariance (plus the requirement of 
2nd order e.o.m. that leads essentially to Lovelock invariants &Co). We refer the curious 
reader to [2, 3] for further reading on these aspects. 

As mentioned, galileon theories are ubiquitous in the literature and have already been 
employed to probe different scenarios: modified gravity, but also early universe cosmology. 
This is certainly the case of the work in [6], where single field galileon theories where shown 
to be compatible with an inflationary phase and to become the leading contributions to the 
inflationary dynamics in a specific, interesting, regime. 

Assuming for a second that there is indeed such a possibility also for the case in question 
here, Eq. (2.10), it is clear why such a theory looks interesting: we are presented with, besides 
the usual kinetic term, a single quartic Lagrangian, one coupling, which presumably can lead 
the dynamics of fluctuations during an inflationary phase. 

If, upon switching on fluctuations, we were to find a sizable region of the parameter 
space where the third order action is suppressed (such an occurrence is clear-cut fine-tuning 
in a scenario with £3^0 and no well justified symmetries), we would have landed on a stable, 
predictive, theory with a naturally large four-point function and a parametrically small three- 
point function, that is /nl- This is certainly interesting in light of the upcoming analysis of 
Planck results (see also [30, 31]). 
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3 Background Analysis and Inflating Solution 

We will show below that our model supports an inflationary phase and that, thanks to the 
inherited SO(N) symmetry, it can be essentially treated as a single field model for the pur- 
poses of this work. In order for us to be able to trust the inflationary theory inherits all the 
desirable properties described above we will proceed in steps [6] . 

- We first consider an Mpi — > oo dS limit and show that there the background is compatible 
with an inflationary solution and that the shift symmetry is broken, so that inflation can 
end. 

- We then go on to show that the potential chosen, although it does break galilean invariance 
as well, it is not renormalized [32-34]. 

We start with: 



One can make immediate use of the SO(N) symmetry by choosing to specialize along the 
direction of one single field in all that follows. To be sure, the consequence of the multifield 
dynamics (i.e. the conversion of isocurvature modes into adiabatic ones) must not be ignored, 
it is also true that they depend on the specific isocurvature-adiabatic conversion mechanism 
and our focus here is not on such matters. 

It is instructive to consider the Mpi — > oo limit obtained by keeping the background fixed to 
de Sitter and by adding a potential of the form V = Vo — A 3 |-7r|. The potential is needed to 
(eventually) break shift symmetry, essentially for a graceful exit, and we chose its form here 
such that it is simple and such that the renormalization properties are safeguarded. Going 
beyond a mass term (i.e. V ~ ir n , n > 2) is not possible for these reasons [6]. 

Being in curved space, we must now of course have a prescription for the covariant 
galileon theory. As it happens a naive d — > V covariantization does not work, as second order 
e.o.m. would be lost. On the other hand, a covariantization procedure has been worked out 
in [19, 20] and that immediately applies to our £4. The action for the background field 7t(t) 
is: 



Following [6] , we now point out several crucial ingredients: 

- the shift symmetry of the potential is not broken in this limit as any piece linear 
in 7f (think of \J—g A 3 n ) is a total derivative. We conclude that any breaking of the shift 
symmetry is Mp\ suppressed. 

- the equation of motion for 7f is compatible with the 7f = const solution. The latter is 
crucial for the non-renormalization properties of the theory: indeed whenever one states that 
galileon theories are not renormalized this is true provided higher order derivatives such as 
7f, Vr are zero or kept small, otherwise other terms would emerge and our analysis could not 
be limited to terms such as the ones in Eq. (2.2) 




(3.1) 




(3.2) 



- 5 - 



Precisely because in this limit the shift symmetry tt — > tt + const is conserved (the 
potential is a total derivative and everything else has derivatives acting) one can derive a 
current and identify the solution for tt from the current conservation equation. Here we shall 
be content with the e.o.m. from the above action, which reads: 



so that tt = const is indeed a solution for tt and therefore the renormalization properties, at 
least in this Mpi — > oo, dS-fixed limit, stay intact. This is not necessarily true later, when 
one accounts for gravity and the galileon symmetry itself is broken, but in that regime the 
breaking will be parametrized by A/Mpi, which we can safely assume to be small. 

In other words, what we have been after in this section can be summed up as show- 
ing that everything works out nicely in this limit and therefore, whatever happens in other 
regimes, will be at the very worst parametrized by K/Mpi (H/Mpi being the other possibil- 
ity) and so should be under control. We have shown that the shift symmetry is preserved 
and argued that considering tt = and the specific choice of potential [6], also the non- 
renormalization theorems still apply. 

Having established that the galileon terms are well behaved, we must find out in what 
regime they are relevant. Indeed, already at the background level, we see that 

r «2- 9 2 # 2 -2 ,o A \ 

where Z = Htt/A 3 . It might not be extremely suggestive in this form, but the less restrictive 
theory without SO(N) symmetry where £3, £5 survive (that is, the explicit expression for 
Eq. 2.2) leaves no doubt that Z is acting like a coupling constant for galileon theories and 
therefore: as long as Z < 1 these theories have little new to say but whenever Z is greater 
than unity we depart from canonical inflation and galileon non-linearities must be accounted 
for, they are indeed the leading contributions. 

It is clear at this stage galileons non-linearities are leading term in the action just above, 
not so much that they might be leading with respect to metric perturbations as well. In order 
to show that we need a quick detour into the effective field theory of inflation (more precisely 
"of fluctuations around an FRW solution") which was set up in [21] and further generalized 
in [22] (for non-Gaussianities produced within this approach see also [24-27]). 

On the E. F.T.I. Approach and Why It Is Safe to Neglect Metric Fluctuations 

In [21] the authors assume they have a stable effective theory around an FRW background 
and work out the most generic form that the perturbations can take. They first choose 
to work in a unitary gauge where the scalar mode 5<p = (Sir for us here) is "eaten by the 
metric" . This choice breaks time-reparametrization invariance and therefore the most generic 
theory will have only space diffeomorphisms as a symmetry. These are in short the premises 
that lead to the following Lagrangian: 



As the free g 00 indices and the presence of the (essentially 3D) extrinsic curvature K^ u signal, 
the action is indeed only invariant under space diffs. The crucial step now is to reintroduce 




(3.3) 




(3.5) 
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full space-time reparametrization employing the so called Stueckelberg trick: one promotes 
the parameter that describes time reparametrization to a field 7r (not to be confused with 
our 7r) which now has a specific gauge transformation and reinstates full diffeomorphisms. 
This all prescription is ideally cast for a decoupling limit: we want the dynamics of tt to 
decouple from gravity. To this aim, it suffices to consider the terms in the action that mix 
5g, 5fr. After the field promotion, one typical term will have the following form 

M^Hg 00 M^H((1 + fr) V° + 2(1 + Z)g 0i diTt + g ij d^d^). (3.6) 

Focusing on the quadratic part of the first term, upon canonically normalizing both (ir c = 
MpiiJ 1 / 2 ^ and g°° = M P ig 00 ), one finds 

M^H(9°° + 2^°° + *V°) = & + 2H 1 ' 2 * c g™ + Hg™. (3.7) 

It is clear that, the 7f-only term having more derivatives, it will lead over any mixed term for 
sufficiently high energy. In this specific example: 

Emix > H 1 / 2 , (3.8) 

where one should note that the value of E m i x can change if the leading kinetic term for fr is 
different as is indeed the case in ghost inflation [37, 38]. This might surely appear to be a 
different approach than the one we take here, but it easy to intuitively spell out the dictionary: 

- if we where to neglect metric perturbations, it is well known that our 5n would be linearly 
related to the curvature fluctuation 

- also, at first order in slow roll, neglecting 5g, the following relation holds [21]: 

C = -Htt, (3.9) 

that is, whenever the metric fluctuations are switched off, the perturbations of the scalar 
degree of freedom that drives inflation is, in both these languages, proportional to C- It is 
then intuitively clear that there must exist for our Sir, as much as it does exist for tt, a regime 
in which metric fluctuations can be safely disregarded. It it in this regime that we choose to 
work below. 

One might well ask at this stage why, if both the reasons and the procedure whereby 
we neglect the metric fluctuations are so clear-cut in the setup of [21] , we chose to work 
in our set up instead. The motivations are to be found in the previous section where the 
background analysis was instrumental in tracking down the crucial properties of the theory 
as we have moved in steps towards the more realistic model which accounts for gravity. In 
order to do it, one needs to have a hold of the specific theory, a theory of its fluctuations is 
often not enough 1 . We return to a related point in the Discussion section. 



1 It must also be said that, with reasonable additional hypotheses, a (possibly stable) and very interesting 
theory of "galileon inflation" in a setup similar to [21] can indeed be constructed, see [7]. 
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4 Non-Gaussianities 



Having written down the theory and discussed its property at length, we move now to calcu- 
late observables such as higher order correlators of the curvature fluctuations £. The current 
bounds on the amplitude and profile of these quantities are, as well known, soon-to-be made 
more stringent by the analysis of the data provided by the Planck mission. 

Before presenting the actual trispectrum calculation, we pause here to note that, despite 
what one might naively expect, a null cubic Lagrangian in the full theory, does not easily lead 
to a small three-point function. We can see this quickly by looking at the full Lagrangian 
and using an estimate of the main contributions to the three(four)-point function. Consider 

" f( V7r ) 2 " g(Vvr) 2 {(nvr) 2 - (V M V*7r)(V»W) - \r(Vtt) 2 } , 

the a4-regulated main bispectrum contribution will be of the type 

S 3 4 ~ ^d5-K(aSTT) 2 ~ ^H 5 5tt 3 (4.1) 
use Sir ~ n£/H, to convert 5ir into the proper observable, the curvature £, 

S^~^H 2 K 3 ( 3 = a A n 2 Z 2 ( 3 - where Z = ^-. (4.2) 

Proceeding in the same way for the trispectrum, 

~ ^d57rd67r(nSir) 2 ~ ~H e 5ir 4 = ^.^ 2 H 2 ( A = a^ 2 Z 2 ( A , (4.3) 

where for this estimate we have assumed all along that the main contribution to correlators 
comes as usual from the around-the-horizon region where n ~ Hit and din ~ H/c s tt. 
Summing up the estimate, we have 

£ 3 (C 3 ) + A(C 4 ) ~ ^ 2 Z 2 (C 3 + C 4 ). (4.4) 

The fact that the third and fourth order fluctuations have the same coefficient clearly means 
that there is no room for a small bispectrum vs large trispectrum finding. Furthermore, 
it also implies that the relation between the corresponding amplitude parameters /jvx and 
ttvl is the standard one: if this is the case, the value of t^l has to be about 2 five order of 
magnitude larger than /jyx in order to be detectable. This is not what happens in the case 
of Eq. (4.4), as the t^l/ Jnl ratio is much smaller, it stays too small for detection even if 
one is willing to allow as small a speed of sound 3 as c s ~ 10~ 2 , which is far smaller than 
what is allowed in our specific setup. 




With a possibly small Jml (compatible with zero) from Planck data in mind, one might 
have hoped that a stable, predictive and very restricted (by symmetries) model of inflation 
with a null cubic interaction in the full theory would be a good candidate for a naturally 

2 This is also a model dependent statement, so ours here is to be considered a rough estimate. 
3 This value is close to current bounds on c s obtained by the corresponding bounds on /jvl for quite generic 
effective field theory models which comprise a c s C 1 dynamics (see also [18]). 
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large four-point function whilst generating a small three-point correlator (/jvx)- This would 
have granted a detectable tnl- 

For it to happen one needs the (naively) additional background quantity 7f which appears 
in (4.3) to be small so that the three point function coefficient can be smaller than the four 
point function one. Upon comparing (4.2), (4.3) we see that the relation 5ir ~ 7f£/iJ evens 
the score of background quantities thus making the two coefficients essentially of the same 
order. We proceed with the full calculation below, returning to the lesson that can be drawn 
from this fact in the Discussion section. 

We now turn to a detailed analysis of fluctuation. These are obtained by expanding 
around an FRW background, to lowest order in slow roll, with ir(t,x) = vf + 6ir(t,x). We 
determine the wave function from the quadratic action: 



So 



where 



71 



0'2 



tH + 27a 4 Z z 7t 



2^2. 



72 



0'2 



2^2. 



7T + Y&Ol^Z 7T 



(4.5) 



Sir = tt£ . (4.6) 



Notice that, short of indulging into serious fine tuning, the value of c s can safely be assumed 
to be in the 1/2 < c s < 1 interval. This will naturally have consequences for the value of the 
bispectrum amplitude, which cannot be very large. 

The quadratic Lagrangian gives the usual dispersion relation and canonical wavefunction (we 
adopt here the Bunch-Davies normalization). To first order in slow roll, we have: 



u(t, k) 



iH 



i/4~7icp? 



(1 + ikc s T)e 



-ikc s 



with C = H£ 



d 3 xu(r, k)e 



tkr 



(4.7) 



Without much ado we proceed to the calculation of the various contributions to Jnl- This 
has been done for the single- field case in [6] and, considering our use of the SO(N) symmetric 
to obtain an effectively single degree of freedom behaviour, we simply and briefly produce 
here a discussion of the shape-function, and plot it. The formalism used is the Schwinger- 
Keyldish[13-15], so called IN-IN, formalism (see the Appendix for further details). 
The two-point function for the observable £ is defined as: 



(C(r, h)C(r, k 2 )) = (2vr)V 3 ) (k x + k 2 ) P(h 



(4.8) 



we also introduce = Pk /c 3 /(27r 2 ), whose value is found from COBE normalization [36]. In 
the two-point function calculations the interactions enter only at the loop level while here we 
focus on tree level diagrams. One therefore moves to the three-point function, where cubic 
interactions enter at tree level. The part of the action cubic in the fluctuations is: 



d xa^a Z 



.#3 

A 3 



18tf£ a - 12£ 2 <f£ 



uHi(di0 2 + mo 2 dk 



(4.9) 



We are after the three-point correlator shape-function, we find it by isolating a momen- 
tum conservation Dirac delta from the three-point function: 



<C(r,fci)C(r, k 2 )C(r, k 3 )) = (2tt) 3 5^ (h + k 2 + k 3 )B(T, h, k 2 , k 3 



(4.10) 
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The amplitude /jvl is denned as: 



NL 



5 B(T,k 1 ,k 2 ,k 3 ) 

6 \ [ P{ki)P{k,2) + perm.] 



(4.11) 



and then, using invariance under the overall size of the |fc|, one plots the function 



1.2 k 2 x B (h h h\ 



(4.12) 



where the first factor has been added to enhance the differences in the shape- functions profile. 
Whenever a theory has many independ coefficients regulating the interactions one obtains 
different sorts of shapefunctions by tuning the different coefficients accordingly. Not so in 
this case, the SO(N) symmetry, no matter the value of c s (which depends also on a 2 , not 
just 0:4), demands an equilateral shapefunction, we plot it in Fig (1). 



B(l jm) 




Figure 1. The shapefunction peaks in the equilateral (k n jk\ = x n = 1) configuration, as expected. 



As for the value of Jnl m this setup, it will depend on both 04 and a 2 or, which is the 
same, «4 and c s . But, considering that the value of c s in most of the parameters space of 
this model resides 4 in the interval [1/2, l] , the value of /nl is easily shown to be \/nl\ < 10 . 

We now proceed to the trispectrum calculation. The quartic contribution 5 to the fluctuations 
Hamiltonian 6 is now needed: 



Oi^ 4 + 2 5n 



8?5lT 



3 u i 



•2 



(diSn) 2 d]5ir (d^vr) 4 (c^vr) 2 {( d j S ^f ~ cW^c^ 



+ 07 



+ 08 



a 



+ 5 5n' 



(didjd-n-)' 2 



(4.13) 



4 Unless one is willing to fine-tune 02,4 to obtain very small values for c s , one can only assume 
012 ~ on , or C12 > «4 , or «2 < 04 . 

5 It is of course necessary to include so called scalar exchange diagrams contributions from £3, but we 
already anticipate these terms signature will not bring any novelty in terms of shape functions. 

6 We aim directly at the Hamiltonian in this case because, as opposed to the cubic case, Ha ^ —£4 and 
therefore one needs to be more careful [35]. 
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where, to first order in slow roll, C,£,c)7r are simply related by £ = — H£ = —(H/tt)57t. The 
expression for T-La with the explicit O n operators it provided in the Appendix. 
The general formula for the calculation of higher order correlators in the IN-IN formalism is 
given by: 

(m 1 Ck 2 c k3 --.CkAm) = (o\f{e i ^ d ^ 

(4.14) 

we spare the reader the details of the analytical results of the calculation in the main text, 
focusing here on the plots of the shape function (also known as form factor); we report them 
in the Appendix. The trispectrum in Fourier space is also defined as: 

(C(r, fcx)C(r, fc 2 )C(r, fc 3 )C(r, h)) = (2^) 9 P ? 3 6& ( £ kA J] ^T(h, k 2 , h, k 4 , k 12 , fe M ) 

\i=i J i= i K i 

(4.15) 

where T is the form factor we will plot. Notice that the four-point correlator shape- function 
depends on six parameters so one has to choose different ^-configurations for a plot. We 
choose them here in order to ease the comparison with existent literature with specific atten- 
tion to the results of other stable (in the sense of effective field theory) inflationary models, 
such as DBI inflation. Our main reference for comparison will be [39]. 
Pictorially, the k-configurations can be understood by looking at the tetrahedron below: 




Figure 2. 

Equilateral configuration: h\ = k% = k% = k 4 , plotting kvz = \k\ + k 2 \,ku = \ki + ki\. 

Folded configuration : ki2 — > 0, k\ = k 2l k% = fct, plotting k^/ki, k 4 jk\ . 

Specialized planar limit: k\ — fc 3 = £44; k 12 = /(&!, &2, £4) , plotting k 2 /ki , k 4 jk\. 

Near double squeezed limit: k 3 = k 4 = fc 12 ; k 2 = g(ki, fc 3 , k 4 , fc 14 , /c 12 ) , k 12 — > 0, plotting k\ 4 , k 4 



What one plots is the T form factor, except in the last configuration, where T • Y\ i is 
preferred, for an enhanced effect. We report below on the plots we obtain for our model. In- 
terestingly, both in the equilateral and the double squeezed configurations the shape- function 
is distinctively different from its P(X, ir) counterpart. The most striking differences are de- 
tailed in the caption of Fig. (3). 

As of today, we do not have as precise a bound on any given trispectrum configuration 
as we do on the bispectrum side (f NL £ { fj^Li ft^Li //vzf 10 }) but, if we were to reach such 
a precision, this shape-functions differences could be translated into something immediately 
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testable. As for the other two configuration, they do not carry as much information as the 
first two in that there the shapes are hardly distinguishable from the P(X, ir) ones. We 
report them in Fig. (4) for completeness. 

All the shapes plotted are obtained for the case of c s ~ l/\/2, in each configuration. 




?..n o.o 



Figure 3. The equilateral configuration on the left strongly differs from the DBI result in that, for 
example the (0, 0) value is finite, non-zero and it is actually the peak of the shape function. This is 
to be compared with a vanishing T at (0, 0) in P(X, ir) theories. 

Also, the double squeezed configuration plot on the right shows a finite non-zero shape function in 
the (fci2 — ^ 0, y = k\i = 1) limit. This is only to be found in contributions to the trispectrum coming 
from third order interaction terms in P(X, it) models. On the contrary, what is plotted here is a so 
called contact interaction contribution, thus coming from the fourth order perturbations. 




Figure 4. The folded configuration form factor is plotted on the left. The planar limit shape is on 
the right side. These configuration do not have particularly distinctive signatures as compared to the 
corresponding plots in [35] 

What we have here then is a quite distinct model of galileon inflation which shares 
with the related co-dimension 1 model [6] all the essential stability features. The additional 
SO(N) symmetry has implemented a theory with only 02, 04 as parameters and resulted in 
a bit less flexibility when it comes to small c s and correspondingly large /jvx- Nevertheless, 
the model accommodates for a /at/, < 10 and has quite intriguing trispectrum signatures. 
We further expand on these point in the conclusions. 

We now offer some comments on the small /nl vs large ttvl possibility; this is a topic 
which, when looked at from a quantum field theory perspective, is interesting beyond the 
cosmological perturbations framework. 
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5 Discussion 



As we have see above, this symmetry dictated galileon model is quite interesting on its own 
as it produces signatures in the trispectrum shape function which clearly distinguish it from 
the P(X, it) (chief among which, DBI inflation) predictions. 

Precisely because of this symmetry, the full cubic Lagrangian is absent from the onset. 
If such an occurrence were to even mildly (one only needs or maybe prefers a parametri- 
cally small three-point function, rather than a null one) propagate itself into the action for 
perturbations, it would be quite intriguing. This is so because the consequently small /jvl 
parameter would imply a t^l at close observational reach. If Planck data will not rule out 
a Inl compatible with zero or if a small /tvl is the most favored one, the very next natural 
step will be to deepen the study of trispectrum shape- functions and tighter bounds on t^l 
would be in sight 7 . 

How then does one obtain in a natural way the combination of a small /jvz, with a 
large tnl^ There are some ideas in the literature, most notably [30, 31] (see [26] for an 
implementation of these in the study of the trispectrum of very generic effective theories), 
which all approach inflation within the effective field theory of inflation methods of [21]. 

Assuming that there exist such a clear-cut symmetry such as e.g. 5tt — > —5tt of [30], 
immediately serves our purpose, but plenty of work is still ahead if one wants to arrive at 
the full effective theory. Indeed, it is necessary to reverse engineer the fluctuations action; 
the process is not easy and certainly not unique. What appears to be a straightforward 
symmetry in the theory of fluctuations, has a very contrived expression in the full one [40]. 

This is a quite generic issue in inflationary setups: 

- there are models which are openly phenomenologically inspired and so are not thought for 
nor suitable to analyze issues such like unitarity and stability of the theory. 

- there are theories which are unitary and stable (e.g. [4, 6]) and therefore are known in their 
full, background independent, form. Their predictions are about to be put to more stringent 
tests by Planck data. 

- there is an effective theory of fluctuations approach [21], which is extremely powerful in 
reading off what happens at the fluctuations level and mapping it onto specific operators 
in the Lagrangian £(5ir n ). In this approach one rightfully assumes that there is an FRW 
background to expand around and cleverly uses it to gain calculational advantage. The clear- 
cut dictionary between non-Gaussianities and specific independent operators is ideally suited 
for the task of predicting almost all possible non-Gaussian scenarios that data might show. 

However, it is quite hard to further map the C(5ir n ) symmetries into their full C(iT n ) 
counterpart. 

- an effective background independent field theory of inflation has also been put forward 
[23]. The galileon model discussed here nicely fits within this setup (DBI inflation goes even 
further in that it can claim UV finite status). It is in this approach that the properties of 
the full theory are more manifest and, also, where the theory is in a form more suitable to 

7 Strictly speaking this is, of course, a model-dependent statement. 
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guess any corresponding UV finite embedding. 

But, as we have seen, it is hard to implement in a simple way seemingly straightforward 
requirements from the non-Gaussianities endpoint of the theory. Further inputs might mate- 
rialize with the upcoming release of Planck data and it might soon be of utmost importance 
to create a dictionary between this approach and the one of [21] if we are after a stable model 
of inflation which is not accommodated by any of the models in e.g. [4, 6]. 

Even if in cosmological setups Jnl in all its probed disguises (local, equilateral, orthog- 
onal etc..) turns out to be accounted for by existing models, this 



dictionary is an interesting quantum field theory issue by its own [40]. Everytime one has, 
as is frequently the case, a window on perturbations around a specific background and a 
symmetry is apparent in that setup, it is not at all obvious how to retrace that into the 
full theory (when available). We plan to expand on this and provide concrete examples of a 
dictionary in [40]. 

6 Conclusions 

We have studied a specific inflationary model in a regime where non-linear galileon interac- 
tions play a leading role. This theory stems from the study of the higher dimensional brane 
constructions origins of galileon interactions, much in the spirit of the original DGP proposal 
[1]. Whenever the dimensionality of the (bulk -brane) space is even, if the galileon interac- 
tions are to preserve their essential flat space and co-dimension 1 properties, an additional 
SO(N) symmetry arises which greatly restricts the allowed building blocks in the theory. 
Indeed, in the AD (3-brane) case, only one coupling regulates the interactions. 

Having a smaller parameters space translates into less flexibility when it comes to the 
value of the speed of sound c s and non-Gaussianities in general. However, we have shown 
that the bispectrum amplitude can go as high as Jnl i$ 10 (possibly more, but we stick to 
10 without resorting to fine tuning) and peaks in the equilateral configuration. 

The trispectrum shape-function is quite intriguing in that it is distinctively different 
than the corresponding P(X, n) results in more than one momenta configuration: equilateral 
and double squeezed. 

We stress here that the galilean origin of this theory endows it with second order equa- 
tions of motion (unitarity at the quantum level) and non-renormalized interactions. The 
latter makes sure the number of leading coefficients in a specific regime stays the same, an 
issue which is typically highly non trivial in generic inflationary models, and thus makes the 
theory properly predictive. These are the grounds on which we have chosen to compare the 
predictions of the model against P(X, n) theories, with especially DB /-inflation in mind. 

Despite in this model all odd (in the scalar field ir) Lagrangian contributions are null 
in the full theory, we have shown that there is no room for a small bispectrum vs large 
trispectrum combination to occur. This fact has prompted some comments on how to realize 
such a scenario and what it means to translate symmetries in the fluctuations Lagrangian 
into their full theory counterparts. This will be essential if Planck data suggests such a 
possibility and represents an interesting question on its own. 

Besides upcoming enhanced sensitivity on cosmological observables, a way to navigate 
the ever-growing space of inflationary models compatible with data is to be more demanding 
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on the effective quantum field theory side, taking an inflationary model seriously enough 
to request quantum stability, if not UV finite origins. Galileon inflation sits nicely in this 
context and we hope to further investigate related models in the future. 
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7 Appendix 

We report below the analytical result of the different contributions to the trispectrum gen- 
erated by the interactions in Eq. (4.13). 

We organize them by operator and therefore have seven different contributions (two of the 
eight different interactions in Eq. (4.13) have identical analytical expression. 



T, 



(i) 



,(3) 



,(4) 



,(5) 



lie! 

H 12 Q 3 
H 1A 4 



H 14 Q 5 

ifc] 3 



WklUikiJ 

16 kfUih J 

3 ^ 

16 klUiki J 

( 3(k-k 4 Y 
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where fa = fa + k 2 + k 3 + k 4 and 
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« 2 I -O^tf 2 + ^ + - 6^4 + ^5 - I 2 ^ _ ID^ + 2_D^\ 

- » 1 C 2 ef cf C | C° y 



+a f 2430D i7 2 Z 2 — 2160Z) 2 Jj ' z2 _ 756D 3 H 2 Z 2 , 252D 4 Z 2 , 324D 5 Z 2 , 2AD 6 HZ 2 _ 20D7H 2 Z 2 , 108-D8Z 

4 V 1 5 <3 5 <3 <3 <5 ^ 

where the -D n represent all the eight different operators which appear in Eq. (4.13). 
The definition of tjvl reads: 

(^component ~> (2tt) 9 5 3 h)^ T NL „ (7.6) 

i 

to be calculated in the regular tetrahedron limit. 
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